Realistic modellization for the growth of a yeast colony. 
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An evolving yeast colony is simulated by means of a cel- 
lular automaton that takes care of many important features 
of the system under study. A complete survey of the prop- 
erties of the colony is done and a set of scaling relations is 
found, together with the analysis of the critical exponents. 
The Family- Vicsek relation between them is verified to a good 
accuracy. 

A mobility parameter is introduced in order to relate the be- 
haviour of the growth exponent with the temperature and we 
propose that the brownian motion is a key feature to describe 
the bacterial and fungi growth process in diluted solutions. 
The chemical and geometrical correlations are found to play a 
relevant role in the modellization, so that the need of realism 
is emphasized. A new phase space built upon the relevant 
'thermodinamical' variables is explored. 

PACS: 87.23-n, 87.18-h 



In the realm of biology the importance of growth of 
bacterial and fungi colonies Q, [2|, p), Q, as well as 
the tissue formation in hystology || or the evolution 
of species ||, are well established field of development. 
When dealing with epidemiology j?J or front propagation 
Q computer simulations inspired to "lattice animals" or 
Cellular Automata |J are widespread implementations 
for these very different systems. 

Furthermore the origins of life are understood on an evo- 
lutive basis, starting from simple molecules arranging 
themselves in more complicated ones. A primitive "first" 
bacterium in a small puddle, in which nutrient substances 
are dissolved or a brewer that is seeding some yeast in 
his malt liquor (the so-called 'wort' ) to produce beer are 
very different examples of the systems we are intended 
to study. 

We think also to the technological importance of yeast 
in the brewing or baking process, or more in general in 
fermentation processes. 

The Eden model |h]] is often cited when the growth pro- 
cess is discussed: actually the shape of the clusters re- 
mind the observed ones, but the scaling behaviour is far 
from the measured one and hence things are much more 
involved. Furthermore to show the capability of com- 
puter simulation to reach a stable state is cited the Con- 
vey's Game of Life [jnj, that could have something in 
common with the evolution of a bacterial colony, but it 
is far from be realistic. These and others early models 
were often based on simple duplication laws, concerned, 
for example, only with the number of neighbouring un- 
occupied sites. 

From all these considerations we are faced with the prob- 



lem of interpreting correctly and modelling simple biolog- 
ical or ecological systems. The main issue of this work is 
to emphasize the fact that is important to include all the 
relevant details to generate a complex enough behaviour. 
It has been suggested || that chemical correlations are 
necessary to explain the very different morphologies and 
evolutive features of growing living systems. We intro- 
duce in our model 'instantaneous' chemical correlations 
represented by the concentration fields of nutrient and 
waste and we explore systematically a new phase dia- 
gram that relates the number of generation of a colony 
to the above cited concentrations. The mobility of the 
cells, that is a function of the temperature, is shown to 
play a relevant role in the evolution of such a system. 
We measure power laws with the purpose to study their 
modifications when the yeast is allowed to wander ran- 
domly. This work starts from that foundations with the 
intent to model properly the system under investigations 
(Sec. II), to generate accurate patterns and to measure 
some important quantity (Sec. II-III) giving estimates of 
scaling exponents. Conclusions on the results are drawn 
in Sec. IV. 



I. THE SYSTEM AND THE MODEL. 

The complex system we are studying needs to be imple- 
mented taking care of all the necessary details. We will 
consider a Cellular Automaton (CA) defined on a square 
lattice (due to the limitations of memory the maximum 
side length is about 800-1000 units) with periodic bound- 
ary conditions. The simulations start all from a single 
seed placed in the center of the grid. The behaviour of 
every cell depends on the eight nearest neighbours: this 
is the Moore type neighbourhood as defined in [Q . 
The typical yeast life cycle is very simple to mimic and we 
examine the simplest behaviour that a cell could keep: it 
eats nutrient (a generic name that means above all sugar, 
with small percentage of other substances) and it elimi- 
nates waste (that in our context is alcohol) . The specific 
chemical composition of these substances is however not 
important here. 

We define two overall initial variables: the concentration 
of nutrient CJ™' d and the concentration of waste C™' te 
(usually set to zero, but one could wish to explore other 
regions in the C f 00 d — C was te phase space) . These quan- 
tities are not properly concentrations but rather amount 
of solute because they are defined as the product of the 
squared size of the grid (the volume) and the percentage 
(a number from to 1). 
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The reason to introduce such variables is simple: in the 
real situation yeast cells cannot tolerate an alcohol con- 
centration higher than 12%. The colony is "killed" by 
poisoning when this value is reached: this is some sort of 
criticality because the system grows, obeying to certain 
probability rules that we will explain below, consuming 
food and producing waste until the refusals are too much 
to let the colony live. For our purposes the initial con- 
centration of food is very large acting as a reservoir, but 
we can imagine situations in which the colony is " killed" 
(or properly inactivated) by starvation. One may think 
to other strategies to implement food and waste like in- 
troducing a food cell that is absorbed when in contact 
with a yeast cell, thus producing one waste cell. This 
is however not necessary because the typical duplication 
time of a colony is of the order of 30-60 minutes, depend- 
ing on various strains of yeast, on the temperature and 
other environmental features. In that amount of time, at 
room temperature and for small enough systems, the so- 
lution has the same concentration in every point because 
a normal yeast cell is an oval of about 5/xm and hence, 
accordingly with what we have stated before, the maxi- 
mum grid could have a size of about 800 ■ bfim = 4mm, 
some millimeters in the biggest case. 
The larger is the area exposed to the solution, the higher 
is the exchange rate both of nutrients and of waste. 
Hence if all the nearest neighbouring cells are occupied 
there is no way to eat: in such condition the yeast sur- 
vives but is in a suspended situation. If we allow the CA 
to move, as explained below, the yeast cell can restart 
to have a normal life in any moment in the future. This 
takes care of the effect of overcrowding in the evolution 
of the system. 

We can imagine that this movement is characterized by 
a simple random walk of one step in one of the eight sur- 
rounding squares (like the King movement in the game of 
chess). A "mobility" parameter T p is related to the av- 
erage brownian motion of the full points in the grid. For 
every generation a cell is chosen at random and moved 
randomly in one of the eight directions (if it is not occu- 
pied by another member of the community). Obviously 
a completely surrounded one cannot move, but it could 
return to walk in any successive step of the colony evo- 
lution. This is repeated [NT p ] times where N is num- 
ber of individuals and T p is the mobility parameter, that 
becomes an integer when we take the floor function as 
indicated with square brackets. In so doing a given cell 
could be picked out more than one time, and there could 
be cells that are not moved. This is intended to por- 
tray a completely random environment and is in some 
way a new type of modcllization for such a specific sys- 
tem. Similar stochastic automata are discussed with very 
great care for example in Ref. ]l3| ], Jl4[ . 
An important point here is the introduction in our CA of 
a probability rule for replication that depends upon all 
the relevant geometric and chemical details included in 
the present discussion. At a certain generation k, that 
corresponds to a certain time tk, the probability for an 



empty cell f of becoming occupied, because of the mitosis 
of a neighbouring occupied cell, is 

Tl em p C fooditk) (^ waste C waste (tk)) 

P\ r > l k) a (~<init rideath ' 
° food ^waste f 

A ' p ' W 

(1. 1) 

where A is an area factor that is proportional to the 
number of empty neighbouring cells n emp , F is a 'food' 
factor proportional to the overall concentration of food, 
normalized to the initial value. It is one at the beginning 
and decreases during the evolution. W is a 'waste' factor 
that is zero at the first cycle and increases proportionally 
to the quantity of alcohol expelled. The functional de- 
pendence of these factors on the concentrations has been 
taken to be linear, but one can think to other ways of 
implementing F and W . From the definitions above we 
see that < p < 1. 

We denote by c the amount of food absorbed by a sin- 
gle cell in every duplication time, that is equal, for the 
sake of simplicity, to the amount of waste produced at 
the same time. Hence we can write the concentrations of 
food and waste at the time tk as follows: 




Ni is the population at the i-th cycle. 
The whole life cycle is iterated thus generating the com- 
plex dynamics. A typical run consists of: 

1. Seeding of the 'first' progenitor. 

2. Iteration of many generations as long as saturation 
is reached. 

3. Measure of saturation values for radius and width. 

In particular for the point 2: every generation means a 
complete life cycle for the colony. The principal rules are 
the following: 

• Metabolism: every yeast cell eats an amount c of 
food thus producing an equal quantity of waste. 

• Mobility: if T ^ then every cell is moved accord- 
ingly to the stochastic rule exposed in the above 
discussion. 

• Duplication: if the probability p(r,tk) is greater 
than a random number, chosen in the range 0—1, 
a cellular subdivision occurs (mitosis). 
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• Measurement: the number of individuals and 
concentrations as well as radius and width are 
recorded. 

Many runs are typically taken into account, and extensive 
simulations are performed with different sets of parame- 
ters, as will be discussed below. 



II. DISCUSSION. 

Our simple automata are suitable to breed a huge va- 
riety of situation, but the gross feature of the evolution 
in essentially the same: after the well-known initial ex- 
ponentially growing phase the duplication becomes hin- 
dered by the cooperation of diverse mechanism, each 
leading to a small rate of birth of new individuals: over- 
crowding, starvation and poisoning act in the same di- 
rection. Note that this concepts are not only effective 
in yeast growing, but a huge class of ecological and bio- 
logical systems is based on the same foundations, first of 
all the human kind, even if the control mechanisms are 
much more intricately involved and hard to analyze. The 
last part of the evolution reaches a saturation value after 
which the growth stops completely because the alcohol 
concentration turns to the critical value. 



A. T p = case, no thermal motion. 

Let us begin with the discussion of the case with no 
brownian motion Tp = 0, when cells are not allowed to 
move on the grid. The shape is very compact and the 
colony develop in a larger number of generation, until 
the poisoning is achieved. A typical evolution leads to 
the behaviour showed in figure 1. 




where d is the distance from the center of mass (xb,yb) 
and Ng is the number of individuals in the surface region. 
The semi-width of the interface between the colony and 
the ambient solution is: 



(2.5) 



These formulas are written in analogy with the corre- 
spondent quantities in a ballistic deposition simulation 
(see, for example, [@). 

In figure 2 we report the width divided by L a as a func- 
tion of the time divided by L z . The collapse is so good 
that one can hardly distinguish the four different symbols 
used for four different values of the grid size. 
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FIG. 2. Data collapse for different values of the grid size: 
100,200,500,800 (circles, squares, diamonds and triangles re- 
spectively) . 

The exponents a, (3 and z are defined measuring the 
scaling relations: 



W S at ~ L " , W sat ~ t 



t ~ L z 
with critical exponents 



(2.6) 
(2.7) 



Tp=0 

FIG. 1. Shape of a colony at the saturation. The grid size 
is 200 and the temperature parameter is set to zero. The 
resulting shape is very compact and no cell could exist without 
a neighbour (or better without a parent). 

We have measured the number of cellular automata 
at every cycle, as well as the concentration of food and 
waste. The radius of gyration of the cluster is defined as: 



R 



N s 



(2.4) 



T P = 



a 0.961 (± 0.005) 
1.055 (± 0.018) 
z 0.909 (± 0.017) 

The errors considered here are only due to the linear fit. 
The relative importance of the statistical error measure- 
ments in the code could be estimated to be of the order 
of 10^ 3 and it has been neglected. They are not shown in 
the loglog plots in figure above. Incidentally, in the T = 
case, we note that the Family- Vicsek scaling relation (see 
for reference |lq|) holds almost perfectly: 



a 

Z = 1 



0.961 
1.055 



= 0.911(±0.016) . 
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With the fitted values for the exponents we have checked 
the data collapse: in the next figure we report the width 
divided by L a as a function of the time divided by L z . 
The collapse is so good that one can hardly distinguish 
the four different symbols in the next figure used for four 
different values of the grid size. 

Our value are in goo d agreement with the few existing 
issues in literature |l5|]lq |. Vicsek gives an estimate for 
the exponent a exp — 0.78 that is different from ours. 
However it must be said that the geometry of the ex- 
periment is rather different: they check the evolution of 
a colony on an agar-agar solid substrate inoculating the 
bacteria with an infected rod. In this way they don't have 
thermal motion of bacteria in appreciable way (a feature 
that we will address below) and the chemical correlations 
are less pronounced being the 'almost instantaneous' di- 
lution of substances prohibited in such an arrangement. 
Furthermore we want to mention another possibility, not 
previously taken in to account: the system is not really 
two-dimensional. Some cell is budded above the colony 
thus forming a 2 or 3 layers system. This, in our opinion, 
could led to a reduction of some critical exponent. 

B. T p case, thermal motion. 

Now we turn to the description of the movements on 
the grid. In Fig. 3 we show that the brownian motion 
affects the width of the surface, as well as the size of the 
colony. We expect that an 'hot' colony feels much more 
the finite size effects, and its scaling exponents might 
change with the mobility parameter. 




Tp=2 

FIG. 3. Shape of a colony at the saturation when the grid 
size is 200 and the mobility parameter is 2. The resulting 
shape has a diffuse width and there are holes in the interior. 
This is the effect of the brownian motion. 

Obviously the surface width is greater than in the case 
of zero temperature and this is due to the introduction of 
brownian motion. Every cell, and in particular a newly 
budded one, could now move, thus allowing the surround- 
ing cells to be in contact with the environment. This im- 
plies that they may eat and reproduce again. 
Let us open a parenthesis on the role of T p , before the 
analysis of the critical exponents. To relate our mobility 
parameter to the real temperature we perform a number 
of simulations not allowing the cell to reproduce. The 



Einstein's relation jl7j for the diffusion could be used to 
relate the square of the diffusion length < x 2 > covered 
by the solitary cell with the interval of time At that is a 
multiple of the time of a life-cycle: 

< x 2 > = 2D At , (2.8) 

where D is the diffusion coefficient expressed in terms 
of the temperature T, the viscosity coefficient of water 
v = 0.001 (SI units), Na v is the Avogadro number, R 
the gas constant and the radius of the particle r: 

D = — ^— . (2.9) 
6tt NAvvr 

Refer again to jlTj, taking care of inserting the Stokes 
formula in the denominator. Hence one obtains: 

m < x 2 > 3-jtNavV < (AL) 2 > 3nN Av v , 

± rsj f r^> 7*° 

At R n g Rt g 

(2 . 10) 

where the second expression is rewritten knowing that 
the distance over the grid < (AL) 2 > = < x 2 > jr 2 
is measured taking the size of the grid r as a unit, and 
that the time elapsed during the simulation is measured 
taking the period of a generation t g as a unit, multiplied 
by the number of generations n g . This is a quite coarse 
procedure. 

Measuring the quantities for the first ratio and giving 
an estimate of the radius of a bacterial cell are the only 
hurdles to overcome. While the code permits to achieve 
easily this ratio, unfortunately the size of a yeast cell 
could vary between 2 [im and 10 [ivn, thus making the 
expectations for T rather rough. It could be said, how- 
ever, that a T p of a few units could led to a reasonable 
temperature. 

The critical exponents are measured in a couple of case: 





Tp=l 


Tp=2 


a 


0.931 (± 0.012) 


0.910 (± 0.010) 




1.096 (± 0.010) 


1.113 (± 0.004) 


z 


0.850 (± 0.016) 


0.818 (± 0.010) 



A general trend is the reduction of the a exponent with 
the temperature, at the same time (5 increases and z de- 
creases. Even if our geometrical set up is very different 
from the already cited experiment of Vicsek and even if 
we are dealing with diluted solutions rather than agar- 
agar plates, the observed reduction indicates that this 
is an essential feature to reproduce the reported experi- 
mental values. The only problem to solve is the delicated 
relation between the mobility of a typical yeast cell and 
our mobility parameter. We have made attempts in this 
direction and it seems that a T p = 8 10 generates a 
good accordance with the cited experimental value. 
We suggest that the value 0.78, that is somewhere in 
the middle between the Eden model and ours, could be 
obtained if one accepts the role of thermal brownian mo- 
tion. The data for a systematically tend towards lower 
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values with increasing mobility parameter. 

Even in these cases the Family- Vicsek relation is verified, 

being the calculated z very similar to the measured one: 



a 



- = 0.849(±0.013) 



T = 2 ^ = 0.818(±0.009) 

The Family- Vicsek scaling relations seem to be verified 
to a good accuracy and even in these cases a data col- 
lapse similar to the previous one has been verified. 



III. PHASE-SPACE WANDERING. 

The chemical correlations established by the diffusion 
of substances in the solution is studied by changing the 
initial parameters. Let us play with a phase space in 
which the relevant axis are N,Cf 00 d and C waste , fixing 
the mobility parameter to zero for semplicity. Two re- 
gions can be identified accordingly to whether the colony 
is able to live or it must die. The localization of the 
surface between the 'life-zone' and the 'death-zone' is in- 
teresting. 

The following figure illustrates the dependence of the 
number of individuals with the amount of food at dis- 
posal varied between — 1 for two different system size. 
As one can imagine N growths if the initial amount of 
food is high, but the mathematical formula could not be 
simple to guess. We believe that it growths like a loga- 
rithm. 



However, as one could expect, the number of individuals 
is much more sensitive on this parameters because this is 
the mechanism that actively stops the growth. 
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FIG. 5. {iV — Cwaste }Section of the phase space. Number 
of individuals vs. critical amount of waste. The data are 
displayed for two different grid sizes as indicated in the figure. 

Having this two sections in mind, one can easily 
reconstruct a surface in the three dimensional space 
{N,C f 00 d, C was te}, even if it's not simple to find the 
other section, namely the one in the {C ' food,C was te\ 
plane, because the code does not allow to fix the number 
of cells that is an output. 

Another interesting graph we want to discuss displays 
the behaviour of the number of generation n g with the 
amount of food. Surprisingly when there is a small quan- 
tity of nutrient in the solution, the number of generation 
is greater. 



N 
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FIG. 4. {iV — Cfood} Section of the phase space. Number of 
individuals vs. initial amount of food. The data are displayed 
for two different grid sizes. 

The same has been done by keeping the amount of food 
fixed (and equal to 1) and changing the critical value for 
the waste from 0.12 to 0.96, for two different grid sizes. 
The result are very similar, but the slope is less steep. 
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FIG. 6. Number of generation that a colony develops before 
being killed by poisoning. The x axis is the initial amount of 
food. The grid size here is 400 and the killing amount of 
alcohol is 12%, as usual. 

This is due to the fact that a few individuals are al- 
lowed to live, but not to reproduce thus slightly consum- 
ing the food storage and producing only a few poison. 
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For intermediate values of food amount there are oscil- 
lations in the number of generation, probably due to the 
effect of balance between the two regimes. 
A similar, but much less impressive plot could be done 
by changing the killing amount of waste and recording 
the number of generations. 
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FIG. 7. Number of generation that a colony develops before 
being killed by poisoning. The x axis is actually the amount 
of alcohol that stops the growth. The grid size here is 400 
and the initial amount of sugar is 100%, as usual. 

In this case the curve is almost linear and very easy 
to discuss: if a yeast cell could bear a larger concentra- 
tion of poison in the solution, the number of generations 
growths linearly. Different yeast strains have different 
tolerance with respect to alcohol concentration and our 
model might be too simple to be of help when very differ- 
ent species are considered. We think that the behaviour 
that we have discussed depends, to a large extent, on the 
functional form of Eq. (1.1). 



IV. CONCLUSIONS. 

We have depicted a typical scenario that holds for a 
large number of complex systems involving living organ- 
ism. An accurate modelling has been achieved by means 
of a cellular automaton defined on the square lattice. The 
generated patterns as well as some important critical ex- 
ponents have been achieved. The role of the brownian 
motion has been discussed and has been outlined the 
importance of a much deeper understanding of the be- 
haviour of the exponents with respect to the temperature 
that seems to influence the dynamic in a cumbersome 
way. 

It is very probable that a proper understanding of the 
movements of cells in bacteria or fungi colonies could be 
the key feature in reproducing the experimentally ob- 
served scaling behaviour. 

We suggest to biologists or physicists interested in this 
field to undertake simple experiments to measure critical 



behaviour of colony growth in solutions of nutrient at dif- 
ferent temperatures with the intent to better understand 
the role of the brownian motion in this context and to 
check our results. 

We remark the fact that this kind of evolution with 
'instantaneous' chemical correlations and thermal move- 
ment is much more akin to the real one: industrial or 
natural fermentation process and evolution of life very 
often occurs in a concentrated solution at room temper- 
ature rather than on an agar plate. A further step will 
be to implement a new geometrical arrangement, similar 
to the Vicsek experiment one, with the goal to reproduce 
the scaling exponents, and to study the possibility of a 
multilayer colony. 

The interest for this kind of studies could be broader 
than what we have depicted: many other features arising 
from different chemical or biochemical mechanisms could 
be used to introduce more complexity and new patterns. 
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